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Abstract. Sonoluminescence is the intriguing phenomenon of strong light flashes 
from tiny bubbles in a liquid. The bubbles are driven by an ultrasonic wave and need to 
be filled with noble gas atoms. Approximating the emitted light by blackbody radiation 
indicates very high temperatures. Although sonoluminescence has been studied 
extensively, the origin of the sudden energy concentration within the bubble collapse 
phase is still controversial. It is hence difficult to further increase the temperature 
inside the bubble for applications like sonochemistry and table top fusion. Here we 
show that the strongly confined nobel gas atoms inside the bubble can be heated very 
rapidly by a weak but highly inhomogeneous electric field as it might occur naturally 
during rapid bubble deformations. An indirect proof of the proposed quantum optical 
heating mechanism would be the detection of the non-thermal emission of photons 
in the optical regime prior to the light fiash. Our model implies that it is possible 
to increase the temperature inside the bubble with the help of appropriately detuned 
laser fields. 
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Figure 1. Time dependence of tlie driving sound pressure and of the bubble radius 
in a typical single-bubble sonoluminescence cycle. Point A marks the beginning of the 
collapse phase in which the bubble approaches its minimum radius very rapidly and 
becomes thermally isolated from the liquid. In point B, the temperature within the 
bubble is significantly increased and a strong light flash occurs. Point C denotes the 
beginning of the expansion phase in which the bubble oscillates around its equilibrium 
radius until it regains its stability. 



1. Introduction 

In 1934, Frenzel and Schultes [1] accidentally discovered a phenomenon that later became 
known as multi-bubble sonoluminescence [21 E]. In order to speed up the development 
process of pictures, they applied ultrasonic waves to a tank with a photographic fluid. 
What they observed were tiny imploding bubbles emitting light at low intensity. Interest 
in this phenomenon increased again in 1989, when Gaitan et al. [1] were able to 
produce the cavitation of a single bubble. Optimising the experimental setup, they 
created a stable bubble whose radius changed periodically in time. In each cycle, 
the bubble suddenly collapses and emits a sharp light pulse. In its new form, the 
phenomenon became known as single-bubble sonoluminescence. Later, the groups of 
Putterman [5l[6l[7l[8l[9l[T0l[IIl[T2], Suslick [H [151 [HI [13 [Tg [19] and others 
[20l [2n [22l [23l [241 [25l [26l [27] perfected these experiments. Luminescence of a cavitating 
bubble has even been induced by pulsed laser excitation [28] . 

Fig. [H shows a typical single-bubble sonoluminescence cycle and indicates the 
relevant time scales [29l [30] . The period of the applied sound wave is around 60 /zs. The 
time dependence of the bubble radius is in good agreement with the laws of classical 
physics and can be described by the Rayleigh-Plesset equations ^31j. For most of the 
cycle, the bubble radius increases isothermal. Each expansion phase is followed by a 
rapid collapse phase. Close to point A, i.e. about 20 ns before the minimum radius is 
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reached, the accelerating bubble wall becomes so fast that the bubble becomes thermally 
isolated from the liquid. Close to its minimum radius of about 0.5 fim, i.e. between 
points B and C, the bubble might be filled with up to 10^ noble gas atoms. At this 
point, a rapid increase of the energy density occurs which is accompanied by the sudden 
emission of light. In case of argon atoms, the light flash lasts for about 40 ps. Afterwards 
a re-expansion phase begins in which the bubble oscillates around its equilibrium radius 
until it regains its stability. 

Detailed measurements of the spectra of the picosecond light flash between 
point B and C in Fig. [1] have been made. Associating the continuum underlying 
these sonoluminescence spectra with blackbody or Bremsstrahlung radiation indicates 
temperatures of at least 10^-10^ K inside the bubble [3 El [IH [121 [32l |33] . It is even 
possible to observe light emission in the ultraviolet regime which hints at temperatures 
of about 10^ K [6]. Noteworthy is the discovery of sharp emission lines in the optical 
regime [T71 [291 [32] • These indicate the population of highly excited energy eigenstates 
of noble gas [IHl [161 IlHl US] and metal atoms [H] which cannot be populated thermally. 
These excitations prove the formation of an opaque plasma core inside the bubble and 
have been observed in several single- and in multi-bubble sonoluminescence experiments 
[m [321 [31] . Recent experiments [151 [161 [HI [19] already verified the existence of a dense 
plasma inside the bubble. 

The light flash at the end of the bubble collapse phase has been attributed to surface 
blackbody radiation [9l Ull [13], neutral and ion Bremsstrahlung [35l |36l ISZl 138] . 
collision- induced emission [391 SO], quantum vacuum radiation [HI |12], and other 
thermal [TOl [lH [H] and non-thermal processes [ISl [IHl [1?]. Other authors describe the 
plasma temperatures inside the bubble by a converging spherical shock wave [35l [48]. 
All of these models reproduce the observed single-bubble sonoluminescence spectra 
qualitatively. Nevertheless, the origin of the sudden energy concentration during the 
final stage of the bubble collapse phase in single-bubble sonoluminescence experiments 
remains a mystery [32l [51] . A valid theoretical model needs to include a mechanism for 
the formation of a plasma, as well as a mechanism which can increase the temperature 
of the plasma even further by at least one order of magnitude. This mechanism needs 
to be able to operate in a solid state-like environment and on the very small length scale 
given by the radius of the bubble. 

In this paper we show that a relatively weak but highly inhomogeneous electric 
field can cause a coupling between the quantised motion and the electronic states of 
strongly confined particles. When combined with spontaneous emission, this coupling 
results in very high heating rates. The origin of this heating lies in the fact that the 
creation of a phonon is more likely than the annihilation of a phonon. This transfers the 
particles to higher and higher temperatures. We then point out that the electric field 
required for this process might occur naturally in sonoluminescence experiments during 
rapid bubble deformations. Moreover, the noble gas atoms inside the bubble form a 
van der Waals gas which reaches its van der Waals core when the bubble approaches its 
minimum radius [291 [22]. They hence experience strong trapping potentials such that 
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their motion becomes quantised. Furthermore we remark that noble gas atoms have a 
stable enough level configuration to undergo a strong heating process. The quantum 
optical heating of nobel gas atoms might therefore play a crucial role in single and 
multibubble sonoluminescence experiments. 

Our model does not contradict current models for the description of sonolumines- 
cence experiments but might explain certain controversial aspects of this phenomenon. 
In addition to pointing out a previously unconsidered heating mechanism, it clarifies the 
role of the noble gas atoms. Indeed, sonoluminescence experiments require a relatively 
high concentration of them inside the bubble [291 ES] . Moreover, our calculations show 
that the inhomogeneous electric field accompanying the heating process would result in 
a non-negligible population of excited atomic states and the non-thermal emission of 
light in the optical regime. The detection of a certain background radiation prior to 
the light fiash during the bubble collapse phase would hence be a first indirect proof of 
our model. This radiation is expected to occur on a nanosecond time scale and should 
result in sharp emission lines at a typical noble gas transition. Its intensity should be 
proportional to the phonon frequency u and to the number of emitting atoms but does 
not depend on their temperature. 

The model we use for the description of the proposed quantum optical heating 
mechanism is similar to the models typically used to describe laser sideband cooling in 
ion trap experiments [IHlEO]. In these experiments a red-detuned laser field is applied. 
Like the above mentioned highly inhomogeneous electric field, it establishes a coupling 
between the quantised motion and the electronic states of trapped atoms. Based on this 
analogy, we conclude this paper with the proposal to enhance the energy concentration 
in sonoluminescence experiments with the help of appropriately 6/ue-detuned laser fields. 
The observation of a certain dependence of the sonoluminescence phenomenon on the 
frequency (and lesser on the intensity) of the applied laser field would support our thesis 
of quantum optical heating. Increasing the temperatures inside the bubble via laser 
driving could assist sonochemistry in the synthesis of a wide range of nanostructured 
materials and in the preparation of biomaterials [52l EHl [M] and might even help to 
facilitate nuclear fusion [55l [56] . 

There are seven sections in this paper. In Section [2l we explain the basic idea 
behind the considered quantum optical heating mechanism. A theoretical description of 
strongly confined atoms inside a highly inhomogeneous electric field is given in Section 
[3] In Section [H we solve the time evolution of a single atom which is typical for all the 
atoms inside the bubble. In Section [5] we discuss the emission of a certain background 
radiation in the optical regime which accompanies the heating process prior to the light 
flash. Section [6] concludes our discussion by pointing out a mechanism to further increase 
the temperature inside the bubble. We finally summarise our results in Section [3 
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2. A quantum optical heating mechanism 

In this paper, we describe a quantum optical heating mechanism which occurs when 
a weak but highly inhomogeneous electric field interacts with strongly confined noble 
gas atoms. Such a field might occur naturally during the sudden contractions of the 
bubble radius in single-bubble sonoluminescence experiments and during rapid bubble 
movements and deformations in multi-bubble sonoluminescence experiments. Polar 
molecules of the liquid dissolve and ionized species might be trapped inside the bubble 
[291 [32l |57] , thereby leading to a non- negligible inhomogeneous charge distribution. This 
charge distribution is further enhanced by the formation of a plasma which might create 
a significant excess of charge density [151 [HI [HI [El [S2l [M] • Although, the electric field 
gradient inside the bubble might be negligible for most of the time, it might become 
very high on the relevant /im length scale which is given by the bubble radius, when the 
bubble approaches its light emission stage. 

When the bubble reaches its minimum radius (cf. point B in Fig. [1] in single- 
bubble sonoluminescence experiments), the mean distance between the noble gas atoms 
becomes so small that interactions between them can be described by a Lennard- Jones 
potential. The atoms experience an equilibrium between repulsive interatomic forces 
due to overlapping orbitals (also referred to as Pauli repulsion) and attractive forces. 
To model the resulting strong confinement of the atoms inside the bubble, we place each 
of them into a trapping potential. This allows us to quantise the atom motion during 
the collapse phase when the bubble reaches its maximum compression. For simplicity 
we restrict ourselves to the assumption of only one phonon mode per atom and denote 
its frequency by u. We expect that the nobel gas atoms have a stable enough level 
configuration to undergo strong heating processes. 

In the next section, we shall see that the gradient of the electric field inside the 
bubble establishes a coupling between the quantised motion and the electronic states of 
each noble gas atom. To model this we denote the ground state of each noble gas atom by 
|0) while |1) is an excited electronic state. The interaction Hamiltonian contains terms 
that correspond to the excitation and to the de-excitation of atoms accompanied by the 
creation and the annihilation of a phonon. Crucial for changing the temperature of the 
atoms is the presence of a large spontaneous decay rate of the excited state |1) which 
keeps the atoms predominantly in their ground state. Although the described processes 
are highly non-resonant, we show below that they can result in a significant change of 
the mean phonon number per atom. Within a few nanoseconds, the temperature inside 
the bubble can increase by many orders of magnitude. 

Before going into a more detailed analysis, let us describe the proposed heating 
mechanism in a more intuitive way. Suppose an atom is initially in its ground state 
and possesses exactly m phonons. We denote this state by |0,m). Fig. [2] shows the 
immediately relevant transitions for this state which we derive explicitly in the next 
section (cf. Eqs. (jH]) and (fTTI) ). Notice that phonons are bosons with their annihilation 
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Figure 2. Level configuration of a single atom-phonon system indicating the 
immediately relevant transitions, if the atom is initially in its ground state |0) and 
possesses exactly m phonons. As in Eq. ([5]), fl denotes the coupling constant for phonon 
conserving transitions to the excited atomic state |1), while A is due to the electric 
field gradient inside the bubble and establishes a coupling between the electronic and 
the motional states of the atom. Moreover, T is the spontaneous photon decay rate of 
level 1. 



and creation operator b and given by 

oo oo 

b = y/m \m — 1) (m| and b^ = y/m + l\m + l){m\ (1) 

m=l m=0 

with [b,b^ = 1. Consequently, a transition from the state |0,m) into |l,m+ 1) occurs 
with a rate proportional to + 1, while the rate for a transition from |0,m) into 
|l,m — 1) scales as ^/rn. When the spontaneous decay rate of the atom is relatively 
large, such a transition is most likely followed by an irreversible and predominantly 
non-radiative transition back into the ground state of the atom. It transfers the atom 
either into its initial state |0,m), into |0,m — 1), or into |0,m + 1). Since the final 
population in the state with m + 1 phonons is higher than the population in the state 
with m — 1 phonons, the net effect of the described excitation and de-excitation process 
is an increase of the mean phonon number per atom, i.e. heating. 

3. Theoretical model 

Let us now have a closer look at the mechanism which creates the coupling between the 
quantised motion and the electronic states of the nobel gas atom shown in Fig. [2j To do 
so we assume the presence of a weak but highly inhomogeneous electric field and derive 
effective rate equations for the time evolution of a single atom-phonon system which is 
typical for the many atoms inside the bubble. For simplicity, we approximate each atom 
by a two-level system. 
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3.1. The interaction Hamiltonian 

Suppose = |1)(0|, = |0)(1| and r is the position of a single atom inside the 
bubble. Within the usual dipole approximation, the Hamiltonian of the atom is given 
by the dipole interaction 

//int = eD-E(r) (2) 

with e being the charge of a single electron, the (real) atomic dipole moment 

D = Doi a- + H.c. , (3) 

and the electric field E(r). For simplicity, we assume that all field components point in 
the same direction of a unit vector k, which allows us to write E(r) as 

E(r) = ^EfcC^^'^-' + c.c. (4) 

k 

with amplitudes E^ and wave vectors k = fck. Moreover, we consider the motion of 
the atom as quantised. For simplicity, we describe its motion in the k-direction by 
only a single phonon mode with frequency z/ and the annihilation operator b in Eq. ([T]). 
Replacing the position operator r — R accordingly, we find [58] 

k-{r-R) = kAx{b + b^), (5) 

where R is the current equilibrium position of the atom, M is its mass, and 

Ax = {h/2Mvf/^ (6) 

would be the width of its ground state wave function when replacing the square-well 
like potential seen by each atom by a harmonic one. An accurate description of the 
motion of the atom would involve a continuous range of phonon frequencies. Here u 
is the dominating frequency of the populated phonon modes. This approximation is 
similar to the approximation made when describing the light emission from a blackbody 
by a single frequency in order to attribute a temperature T. 

Eq. ([6]) can be used to get a very rough estimate of the typical phonon frequency 
V in sonoluminescence experiments. Suppose, the bubble has a radius of 500 nm and 
is filled with about 10^ argon atoms. Although this might vary by at least one order 
of magnitude, we assume that the volume of each atom is given by Ax^. Comparing 
this volume with V/N where V is the volume of the bubble, we obtain a typical Ax 
of about 4nm. Using Eq. ([6]) and taking as an example the atomic mass of argon into 
account, we find that the phonon frequency v is of the order of 50 MHz. This frequency 
is roughly of the same order of magnitude (or higher) as the phonon frequency in typical 
ion trap experiments [HI [50] . 

Since the atom is well localized within the wavelength of its trapping potential, 
i.e. kAx <C 1, we can now apply the Lamb-Dicke approximation which is routinely used 
in the theoretical modelling of ion trap experiments. It allows us to assume [58] 

exp(ifck- (r-R)) = l + iA;Ax(6 + 6^). (7) 
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Substituting this into Eq. ([2]), we obtain the interaction Hamiltonian 

i7i,t = fi^ {cr- + a+) + M (6 + b^) (a" + a+) (8) 

with the (real and positive) couphng constants 

Q = {2e/h) J2 Doi ■ Re (Eke'''^'^^ , 

k 

A = - {2eAx/h) ^ Doi ■ Im (^E^e^'^'^-^) . (9) 

k 

This Hamiltonian is essentially a Jaynes-Cummings Hamiltonian [SS]- A closer look at 
Eq. (Q shows that A is proportional to the gradient of Q in the direction of the quantised 
motion of the atom, i.e. 

A = Ax k ■ Vr](R) . (10) 

A strong atom-phonon coupling therefore does not necessarily require the presence of 
a strong electric field. It only requires a highly inhomogeneous electric field inside the 
bubble. 

3.2. The master equation 

Dissipation in the form of spontaneous photon emission from the excited electronic state 
|1) is as usual taken into account by the master equation [59] 



p= — -[i^int + ^0 + hub'^b , p] + r[a pa^ — \a^a p 

-\pa+a-]. (11) 

Here hujQ is the energy of the excited state |1) with the spontaneous decay rate V. Instead 
of integrating the master equation, we now use it to obtain a closed set of effective rate 
equations which can be solved more easily. Notice that some of the major quantities in 
these equations are coherences and not populations. 

3.3. Effective rate equations 

In this paper we are especially interested in the mean phonon number per atom and the 
population in the excited state |1). These are given by 

m = {b^b) and X3 = {a+a' - a'a^) . (12) 

As we shall see below, a crucial role in their time evolution is played by the expectation 
values 

Xi = (a+ + a-) , X2 = - a+) , 
Fi = (6 + 6t), Y2 = i{b-b^), 

Y, = {b^ + b^^), Y^ = i{b'-b^^) (13) 
and the atom-phonon coherences 

Z, = ((a+ + a-)ib + b^)) , Z2 = i((a- - a+)(6 + b^)) , 

Z3 = i((a+ + a-)(6-6t)), z,^-{{a--a+){b-b^)). (14) 
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To obtain a closed and relatively simple system of rate equations, we assume in the 
following that 

Uq ^ u, r, Q, A and m ^ 1 . (15) 

Moreover, we approximate the expectation value of operators of the form (Ba^) by 
(B ){o'3), as it applies when the expectation value of B is about the same for an atom in 
|0) and for an atom in |1). Doing so, we find that the time derivatives of the variables 
in Eqs. (IT2l)- (|T^ are given by 

m = AZ3 , X3 = 2{nX2 + AZ2) - r (X3 + 1) , 
Yi = -uY2, Fa = 2AXi + z/Fi , 

Y3 = -2{uY^ + AZ3), n = 2(z/r3 + AZi) , (16) 

while 

X^ = -uJoX2, X2 = -2in + AY^)X3 + uJoX^, 

= -UJ0Z2 , Z2 = -2{nYi + AF3 + 2Am)X3 + cjq^i , 
Z3 = 2A-cJo^4, Z^ = -2{nY2 + AY4)X^ + UJ0Z3 . (17) 

These rate equations hold up to first order in 1/ujq. 

Taking into account that ujq is much larger than all other system parameters 
(cf. Eq. (fTSll ). we find that the rate equations in Eq. ( JTTl) can be eliminated adiabatically 
from the full set of rate equations. Doing so, Xi, X2, and the Z coherences become 
2{n + AY,)Xs 

Ai = , A2 — U, 

2(fiFi + AF3 + 2Am)X3 ^ ^ 
Zi — , Z2 — u , 

in first order in I/ujq. These quasi-stationary state values are reached on the very fast 
time scale given by I/ujq which is typically of the order of femtoseconds. Substituting 
them into Eq. (fT6l) we finally obtain the effective rate equations 

x, = -r{x, + i) , 

Y,^-.Y,, n.iMl±^ + .r., (19) 

LJQ 

Moreover we have 

2A(ny2 + AF4)X3 



UJo 

Notice that the three differential equations (fTOjl decouple from the rest and can be solved 
separately. 
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4. The dynamics of the system 

In the following, we consider two cases. First we discuss the case of time-independent 
system parameters in the strong atom-phonon coupling regime with 

A > and 4A^ > uuo (21) 

and show that the effective rate equations (fT9l) and ( l20l) can indeed result in a very strong 
and approximately exponential heating process. We also identify the most important 
quantities in the time evolution of the system. In actual sonoluminescence experiments, 
the system parameters can change rapidly in time, especially during the final stage of 
the bubble collapse phase. As we shall see below, rapidly changing parameters might 
enhance the heating of the atom inside the bubble even further. 



4.1. Time-independent system parameters 

In the beginning of each sonoluminescence cycle, the atoms experience neither a strong 
trapping potential nor an inhomogeneous electric field. We can therefore assume that 
the coherences defined in Eq. f|T3|) are initially zero and that the atom is in its ground 
state. Due to having previously been in a hot environment, the mean phonon number 
per atom, mg, is in general already relatively high when the trap is formed during 
the bubble collapse phase. Solving Eq. ( fT9l) for these initial conditions and the strong 
coupling regime assumed in Eq. fl2T|) yields 

X3(t) = - 1 (22) 

up to first order in 1/ujq. Moreover, we find in good agreement with a numerical solution 
of the full rate equations (cf. Fig. [3](a)) that 

n(t) = ^[cosh(At)-l], 

Y^it) = - 1^ sinh(At) (23) 

with the (real) frequency A defined as 

4A2 X V2 



A = ^ 1 . (24) 

\UUJo J 

Eq. ([22D shows that the atom remains predominantly in its ground state. Its role is to 
act as a catalyser for the heating process described in Section [2l For times t of the order 
of 1/A, Yi and Y2 are only of the order of 1/uJo- We can hence safely assume that 

nY2 < AY^ and Q < 2Am (25) 



on a timescale of a few nanoseconds. Taking this and Eq. fl22|) into account, we find 
that the expressions for Zi and in Eq. f|T8|) simplify to 



Zi = -2A(2m + l3)/^o and Z3 = -2AY4/UJ0 (26) 



Sonoluminescence and quantum optical heating 11 




1 2 3 4 5 

time t in units of ns 

Figure 3. The Y coherences and the mean phonon number m per atom as a function 
of time for = 10 MHz while Q = 10^ Hz, A = lO^^ Hz, T = lO^^ Hz [BD], and 
Wo = 10^^ Hz obtained from a numerical solution of the full rate equations (fT6|) and 
([T7)) . Good agreement is found with the analytical solutions in Eqs. (|23p . (|^^ . and 
([30]) which are represented by the shaded areas. 
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3 6 9 12 15 



time t in units of ns 

Figure 4. Logarithmic plot of the heating rate rh/m as a function of time in units of 
2A for 1/ ^ 10 MHz and v = 100 MHz and for luq, O, A, and T as in Fig. [31 Again, good 
agreement is found between the numerical solution of the full rate equations which can 
be found in Eqs. (HH) and HT]) and the analytical solution (cf Eqs. ([^ - ([31]) ). 



in first order in l/uo- The remaining variables m, Y3, and Y4 in Eq. (fT6|) hence evolve 
according to the differential equation 

V = Mv (27) 

with V = (m, Ys, Y^)^ and 



/ 



M = — 



A2 





-(2A2 - uujo) 

y 4A2 2A2 - uuq 

Solving this equation for m(0) = uiq and 13(0) = 14(0) = 0, we find 

8A2(2A2 



(28) 



X2, ,2 



mo sinh (At) 



4A2 

Yi{t) = - - — mo sinh(2At) 



(29) 



which is in good agreement with the numerical solution of the full rate equations in 
Figs. [21(b) and[21^c). The mean number of phonons per atom evolves according to 

SA^ 



m{t) = mo + 



mo sinh (At) . 



(30) 



Fig.[3](c) shows that this equation indeed describes an approximately exponential heating 
process. However, as shown in Fig. [H the relative heating rate rh/m is not constant in 
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time but has a relatively strong time dependence itself. At t = 0, we have m/m = 0. It 
then increases rapidly and assumes values well above 2A after a very short time. From 
Eqs. (EZD and (EH]) we see that 

*=-?^ll. (31) 

m Uq m 

For times t of the order of 1/A, the relative heating rate m/m hence approaches 2A, 
since the expressions 2sinh^(At) and sinh(2A)f:) in Eq. (1291) are of about the same size in 
this case. 



4-2. Time- dependent system parameters 

Up to now we considered only time-independent system parameters. However, in actual 
sonoluminescence experiments, system parameters like the phonon frequency u and the 
atom-phonon coupling constant A might change very rapidly in time. The reason for 
this is that A originates from a highly inhomogeneous electric field which depends, 
for example, strongly on the bubble radius. The phonon frequency u also relates to 
the radius of the bubble since it depends on the trapping potential experienced by each 
atom. To take this into account, we now analyse the effect of time- dependent parameters 
z/(t) and A(t). The concrete size of F and fl is less crucial for the time evolution of the 
mean phonon number m, as we can see from the absence of these two parameters in 
Eq. (EHD. 

Let us assume that u and A change on a time scale t which is long compared to 
the very short time scale given by I/ujq. Then the above adiabatic elimination of the 
differential equations f|T71) remains valid and the dynamics of the atom-phonon system 
is given by Eqs. ( fT9l) and ( !20l) . Moreover, we have seen in the previous section that 
Yi and Y2 effectively do not contribute to the time evolution of m, Y3, and Y4. Their 
dynamics is hence to a very good approximation given by Eq. (1271) . In the following, 
we solve this differential equation for rapidly changing system parameters u and A for 
different situations analytically and numerically. 

To do so, we first notice that the eigenvalues Aj of M are given by 

Ai = and Aa.s = ±2A (32) 
with A as in Eq. fl24l) . The corresponding eigenvectors Uj of M are given by 

ui = ^ / 2(z/2 + A2) 



1 



U2,3 




(33) 



V^(z/2 + A2) 

These eigenvectors are in general not pairwise orthogonal. Our calculation below 
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therefore involves also their dual vectors u* which equal 

/ 2(A2 - z/2) \ 



u 



4v^(Az/)2 



v 



Z/2 + A2 





u 



2,3 



16 (Az/)2 



\ 



2z/A 



A2) \ 



(34) 



with u-' ■ Uj = 6ij. This means, the dual vector u-' is orthogonal to the two vectors Uj 
with i ^ j and its overlap with Uj is unity. 

In the following, we write the vector v = (m, Y3, 14)"'" in Eq. ( 1271) as a superposition 
of the eigenvectors of M. Denoting the respective coefficients by q, i.e. 

3 

V = ^ Q Ui , (35) 

i=l 

using the fact that Muj = Aj u,, and taking the time dependence of the eigenvectors of 
M into account, we can then show that 

3 3 

V = ^ Q Uj + Q Ui = ^ CiXi Ui . (36) 

1=1 i=l 

Multiplying this equation with the u-' finally yields 

3 



(37) 



1=1 



This shows that the coefficients Cj evolve according to 



Cl 



C2 



C3 



A^ 




A4 




A2 




A2 


+ Z/2 


A2 




A2 


+ 1/2 



Cl 



C2 



2vVTA4 

5 Z/2 + A2 

5 Z/2 + A2 



IC2 - C3J 



C3 



8 Vz/^ + 
'5 Z/2 + A2 



Cl 



dt 



A 



In 



:Cl 



— In - 

dt 

^Inl- 



2Ac2 



2AC3. 



These differential equations take the effect of the rapidly changing parameters A (which 
depends on u and A) and i' into account and can be used to analyse the dynamics of 
V (i.e. m, 13, and Y4) for times t of the order of 1/A. Unfortunately, their analytical 
solution is in general not straightforward. 

To illustrate that rapidly changing system parameters generally do not suppress the 
above described heating process, we now consider the following scenario as an example. 
Let us assume that A and u change such that their ratio rj, 

v = X/iy, (39) 
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remains about the same. Using Eqs. (jSD, (ITU]) and (^^, we see that rj equals 



V 



2h [k-vn{R)\ 



1 1/2 



(40) 



This means, the ratio rj is constant when the gradient of the electric field VrO in the 
direction k of the quantized motion of the atom experiences the same time dependence as 
the phonon frequency u. Now the differential equations (l38ll can be solved analytically. 
Given the same initial conditions as in Section HTTj i.e. m(0) = mo and ^3(0) = ^4(0) = 0, 
this yields 



m(t) = mo + mo 



[l + V ' 



sinh 



[ drA(r) . 
Jo 



(41) 



In case of time-independent system parameters, this equation simplifies to Eq. (130|) . 
More importantly, the time derivative m of Eq. (HT!) is always positive, independent of 
how A and u change in time. 

It is even possible to show that rapid changes of the system parameters can 
significantly enhance the quantum optical heating process which we analyse in this 
paper. To do so we now consider sudden variations of t] and ask how they affect the 
time evolution of m. More concretely, we assume that the time dependence of rj can be 
described by an equation of the form 

7]it)=r]oe^^'-'"^ (42) 

for times t > to- When combining Eqs. (l24l) . (127|) and (!35l) . we find that 

^ (C2 + C3). (43) 



m = — 



V5 



Let us now have a closer look at the differential equations for C2 and C3. For sufficiently 
rapid parameter changes, i.e. for \k\ ^ A, we can safely assume that r] changes on a 
much shorter time scale than 1/A. This allows us to neglect the terms 2Ac2 and — 2AC3 
in Eq. (138|) . Moreover we assume that either X ^ u ot u ^ X. Then Eq. ( l38l) can be 
used to obtain 



^(C2 + C3) = ±^[lnr7(t)](c2 + C3) 



(44) 



where the plus sign applies when X ^ u and the minus sign when u ^ X. This equation 
can be solved analytically for ?7's as in Eq. fH2!) . Doing so and substituting its solution 
into Eq. (H3ll we obtain 

X{t) 



m{t) 



A(to) 



m(to 



(45) 



From this equation we see that the derivative m remains always positive as long as m(to) 
is positive. For k, > and A <^ z/, for example, the increase of the mean phonon number 
m slows down when the variation of the system parameters occurs. However, it is also 
possible that the rate m increases exponentially in time. This applies, for example, 
when K > and X ^ u. 
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Figure 5. Time evolution of the mean phonon number m in units of the initial phonon 
number ttlq for different sudden variations of A and the phonon frequency v obtained 
from a numerical solution of the differential equations ((38|) . In (a) we have A = 90 MHz, 
v = 1 MHz and both remain constant in time. In (b)-(d) however, A and v change 
suddenly in time as shown with \{t) — 9Qi^{t). Here r]{t) is as in Eq. (|46|) with (b) 
K = 0.9 GHz and At = Ins, (c) k = 5 GHz and At = 0.15 ns, (d) k = -0.5 GHz and 
At = 1 ns, and (e) k = —3 GHz and A< = 0.1 ns. 



The above analytical result is in good agreement with the numerical solutions of the 
rate equations (!38|) shown in Fig. [51 In this figure, we consider four different examples of 
the scenario where r] undergoes a rapid variation within a small time interval {to, to+At). 
In analogy to Eq. ( H2l) we assume 

Vit) = [e^^'-'"^ - 1] r^o^to,to+|A*(t) + [e-'^^*-*-^*) - l] mOt,^At,to+Atit) 

(46) 
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for times t between to and to + At. Here 6t^ t,^{t) = 1 if ti < t < t2 and 6t^ t^{t) = 
for any other t. As expected, Figs. [5](b) and (c) show an enhanced growth of the mean 
phonon number m compared to the case of constant system parameters which is shown 
in Fig. [5](a). In Figs. [5]^d) and (e), the growth of the phonon number m slows down. 
Most importantly, the calculations in this section confirm that rapidly changing system 
parameters can yield much larger temperature changes inside the bubble than the ones 
already predicted in the previous section for time-independent parameters. 



5. Background emission in the optical regime 

For simplicity, we now consider again time-independent system parameters. In this 
case, X3 = —1 to a very good approximation (cf. Eq. (l22l) ) and the atom remains 
predominantly in its ground state. However, more detailed calculations, going up to 
second order in I/ujq, show that a small amount of population builds up in the excited 
state |1) throughout the heating process. We denote this population in the following by 
Pi and notice that 

Pi = (X3 + l)/2. (47) 

To calculate Pi for times t of the order of 1/A, we have a closer look at the differential 
equation for X3 in Eq. ( fT6|) . Solving it via an adiabatic evolution based on the 
assumption of a sufficiently large spontaneous decay rate F, i.e. 

F > z/, fi, A , (48) 

and combining the result with Eq. (1471) we find 

Pi = ^ . (49) 

To calculate X2 and Z2 up to second order in I/cuq? "we now consider the second 
derivatives of the fast evolving variables Xi, X2 and the Z coherences and notice that 
they evolve much faster than all other second derivatives. In fact, they are of the order of 
cUq. This allows us to adiabatically eliminate Xi, X2 and the Z's from the time evolution 
of the system. Using Eqs. (fT7|) and (!22l) . this yields X2 = and Z2 = —^uAY^/uq up 
to second order in I/ujq. Combining these expressions with Eqs. ( |27l) and (l49l) . one can 
then show that 

Pi = -?^m (50) 

to a very good approximation. The signature even of a small but non-negligible 
population in the excited states of the atoms is the spontaneous emission of light in 
the optical regime at a measurable photon rate /. The reason for this is that I is the 
sum of the contributions of all the atoms involved in the proposed heating process. If 
we denote the relevant atom number by N, then I = NT Pi which results together with 
Eq. (I50D in 

/ = ?^A. (51) 
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Figure 6. Intensity of tiie emitted light per atom, I/N, as a function of time for an 



initial phonon number of ttiq = 10 



10 MHz and v = 100 MHz, and u>o, T, A, and 



O as in Fig. [31 Good agreement is found between the numerical solution of the full 
rate equations and fTT]) and the analytical expression in Eq. (|5ip . 



This means, the intensity of the emitted hght I via the 0-1 transition of the nobel gas 
atom is proportional to the phonon frequency u and to the heating rate m but does not 
depend on the mean phonon number m. 

Fig.[6]shows a very good agreement of this analytical result with the photon emission 
rate I obtained from a numerical solution of the rate equations f|T6|) and f|T7|) . Moreover, 
the figure shows that / increases rapidly during the collapse phase of the bubble. We 
see an approximately exponential growth of / as a function of time. Unfortunately, 
the actual emission of a photon might transfer an atom into a state |0') which differs 
from the state |0) of Fig. [2] for example by the spin quantum number. Eventually, this 
results in a complete depletion of the electronic states |0) and |1). Then the coherences 
in Eq. ( fT2l) become zero, thereby ending the heating process and the emission of light 
from the corresponding nobel gas transition. 

Notice that the light emission described in this subsection is not a description of 
the strong light flash at the end of the bubble collapse phase. This light flash and its 
spectrum are dominated by thermal radiation which can be attributed to very high 
temperatures inside the bubble. Here we are interested in a mechanism which might 
contribute substantially to the creation of these very high temperatures. The photons 
which we describe in this section accompany the proposed quantum optical mechanism 
as background radiation. Their detection would constitute a first indirect proof of our 
model. The background photons should be present prior to the light flash during the 



Sonoluminescence and quantum optical heating 



19 



bubble collapse phase. They are expected to occur on a nanosecond time scale, are 
highly non-thermal, and should result in sharp emission lines at a frequency which is 
typical for the noble gas atoms inside the bubble. 

For the experimental parameters in Fig. [6] one can easily estimate that a single 
atom inside the bubble emits between 10~^ and 10~^ photons per nanosecond. This 
means, a single atom emits on average between 10~^ and 10~^ photons prior to the 
picosecond light flash which ends the bubble collapse phase. The intensity of the emitted 
background photons depends strongly on the total number of atoms involved in the 
heating mechanism. Although we expect this contribution to remain small compared 
to the 10^-10^ photons emitted in the optical regime during the picosecond light flash 
[6H [55] , it should be feasible to detect them. 

6. Possible optical enhancement of the energy concentration 

One way to enhance the energy concentration in sonoluminescence experiments is 
to create an auxiliary highly inhomogeneous electric field inside the bubble thereby 
increasing the coupling between the quantised motion and the electronic states of the 
atoms. However, the atoms experience a very strong trapping potential, especially 
during the collapse phase in single-bubble sonoluminescence experiments, as one can 
see for example from the robustness of the spherical symmetry of the bubble during 
each cycle |29l [32]. An analogy between ion trap and sonoluminescence experiments 
therefore suggests another, potentially more efficient approach: namely the application 
of a blue-detuned laser field which should excite the 0-1 transition of the noble gas atoms 
(or of another atomic species with sharp optical transition lines which have already been 
observed experimentally [17]) inside the bubble. 

In ion trap experiments, an ion experiences a strong trapping potential such that 
its motion becomes quantised P9| [50] . Moreover, a laser field is applied and creates a 
coupling between the quantised motion and the electronic states of the ion. In case of 
a red-detuned laser, transitions where the ion becomes excited and loses a phonon are 
far more likely than transitions where the atom becomes excited and gains a phonon. 
Followed by the return of the ion to its ground state via the spontaneous emission of a 
photon, the overall effect is a permanent reduction of the mean phonon number in the 
ions motion, i.e. cooling. 

The similarity of the level scheme in Fig. [2] and the atomic level configuration 
in ion trap experiments [IHl [50] suggests that a laser field applied to the nobel gas 
atoms trapped inside the bubble in sonoluminescence experiments either enhances or 
suppresses the heating process, thereby testing our hypothesis of the presence of tightly 
trapped atoms during the bubble collapse phase. A maximum enhancement of the 
heating process should occur when the laser frequency equals the sum of the transition 
frequency uq and the phonon frequency z/, as illustrated in Fig. [71 In this case, transitions 
between |0,m) and |l,m + 1) are resonantly driven while transitions between |0,m) 
and |l,m — 1) are highly detuned. Additional laser fields should be used to avoid the 
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Figure 7. Possible enliancement of tlie energy concentration in sonolumiscence 
experiments via laser heating. Shown are the relevant transitions in a single atom- 
phonon system, when the atom is in its ground state |0), possesses exactly m phonons 
and a laser field excites the blue sideband resonantly. 



premature depletion of |0) and |1) due to the spontaneous emission of photons into 
states outside the relevant two-level system. 

Of course, sonoluminescence experiments are very different from ion trap 
experiments. There are for example extreme differences in the temperature and in 
the time-dependence of the trapping potential observed by each atom in its finite size 
environment. Ion trap experiments take place in a far more controlled environment than 
sonoluminescence experiments, i.e. at temperatures close to absolute zero. Nevertheless, 
we expect that a laser is capable of driving the transitions in the atom-phonon systems 
shown in Figs. [2] and [71 Possible transitions, which can be used for this purpose, are 
indicated by the emission of light prior to the light flash at certain optical frequencies, as 
discussed in Section |5l We expect that the sonoluminescence phenomenon depends not 
only on the intensity of an applied laser field but also shows a very strong dependence 
on the laser frequency. 

7. Conclusions 

In this paper we discuss a quantum optical heating mechanism which might contribute 
substantially to the sudden energy concentration in sonoluminescence experiments. Our 
model is based on two assumptions. First, we assume a very strong confinement of the 
noble gas atoms inside the bubble during the light emission phase such that their motion 
becomes quantised just before the maximum compression of the bubble in single-bubble 
sonoluminescence experiments and during rapid bubble deformations in multi-bubble 
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sonoluminescence experiments. Secondly, our model requires the presence of a weak 
but highly inhomogeneous electric field. The electric field establishes a strong coupling 
between the quantised motion and the electronic states of each atom. Heating is due 
to the fact that processes which create a phonon are more likely than processes which 
annihilate a phonon. 

Our model does not contradict current models for the description of sonolumines- 
cence experiments, but explains certain controversial aspects of this phenomenon. We 
show that the phonon energy in the bubble can easily increase by a factor of ten or more 
(cf. Fig.[3]^c)). The relation mhu = k^T can be used to estimate the corresponding tem- 
perature increase. In fact, our model can easily predict temperatures well above 10^ K 
inside the bubble. Rapid changes of the system parameters might further enhance the 
proposed quantum optical heating process. Moreover, our calculations emphasize the 
role of the noble gas atoms inside the bubbles. Although predominantly in their ground 
state, they act like catalysers and facilitate the creation of phonons at a very high rate. 

Our model is based on a quantum optical approach which is routinely used to 
describe the laser coohng of tightly trapped ions. To test our hypothesis of a strongly 
confined noble gas atoms inside the bubble during the bubble collapse phase, we propose 
to manipulate them with the help of an external laser field. Such a laser would create 
a coupling between the quantised motion and the electronic states of the atoms similar 
to the coupling created by a weak but highly inhomogeneous electric field. Depending 
on its laser frequency with respect to the relevant noble gas transition (or with respect 
to a sharp optical transition of another atomic species it is therefore expected to 

either enhance or inhibit the sonoluminescence phenomenon. 

The possibility to increase the temperature inside the bubble in typical 
sonoluminescence experiments with the help of appropriate laser fields would be a direct 
proof of our model. A first indirect proof could be the observation of a background 
radiation prior to the light fiash which accompanies the proposed quantum optical 
heating process. Its intensity should not depend on the temperature inside the bubble 
but is proportional to the heating rate rh, the phonon frequency u and the number of 
nobel gas atoms involved in the heating process. The frequency of the emitted photons 
should correspond to typical nobel gas transitions in the optical regime. 
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